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I. Introduction 

Despite some scattered signs of improvement, the world economic situation is still facing a serious jobs 
challenge, with global unemployment above its pre-crisis level and unemployment in the euro area rising rapidly. 
Most developed economies are still struggling to overcome four major issues: the continued deleveraging by 
banks, firms and households holding back normal credit flows and consumer and investment demand, high 
unemployment which is both cause and effect of the lack of economic recovery, third, the fiscal austerity 
responses to deal with rising public debts and fourth, bank exposure to sovereign debts, the financial sector 
fragility (World Economic Situation and Prospects, 2012). Due to worsening economic conditions and drastic 
spending cuts, unemployment in the euro area rises to record levels and it is expected to stabilize and gradually 
decline as growth recovers (The Current Economic Situation - Euro Challenge, March 2012). A long standing 
challenge in macroeconomics is developing a unique and comprehensive model which could explain all relevant 
observations to the level and evolution of unemployment. Since the early 1960, the interest in dynamical 
characteristics of unemployment has been rapidly growing. This study aims to develop dynamic models in order 
to explain a set of facts regarding job flows, unemployment and output. The paper is organized as follows: after a 
brief introduction, the second section brings evidence on the work of Yuri Kuznetsov and Ronald Shone regarding 
the dynamics of unemployment (and employment). The core of dynamics is bifurcation theory which is 
fundamental to systems theory. 

Section 3 introduces basic definitions and theorems of bifurcation theory employed in this paper, Hopf 
bifurcation being the most commonly seen type of bifurcation since the existence of its boundary is accompanied 
by regular oscillations when the parameters are within a neighborhood of the boundary, and where the oscillations 
may damp to a stable steady state or may never damp depending upon the side of the bifurcation boundary on 
which the parameter might lie (William A. Barnett et all., 2008). As within the development of the theory, the 
matching function has been widely used to help model the inefficiencies in labor markets, further on, the study 
extends the framework along with the findings of its most proeminent theorists, Diamond, Mortensen and 
Pissarides. 

Section 4 extends Ronald Shone's approach to the dynamics of unemployment by providing two 
different settings for the matching rate m(u,v). Propositions establishing the existence of Hopf bifurcation are also 
formulated and proved for each case. This section also lays out their phase diagrams and their economic 
application. Central results in this research are the theorems on the existence and location of Hopf bifurcation 
boundaries in each of the considered cases. The main results of the research are summarized in section 5. 

II. The Basic Model 

In order to study the dynamics of labour market, this research is based on the model developed by Ronald 
Shone (2002), who considers that at the ruling wage there is full employment in the sense that the number of jobs 
is matched by the number of households seeking employment. The working population population, N, is 
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considered to be fixed and the number of jobs available to be constant. At any instant of time a fraction S of 
individuals become unemployed and search over firms to find a suitable job. Let / denote the probability of 
finding a job, i.e., the fraction finding a job. At any moment of time, if U is the fraction of the participating 
labour force unemployed, then 

5(1 — u)N — individuals entering the unemploymet pool 
fuN = individuals exiting the unemployment pool 

The change in the unemployment pool, uN, is therefore given by the differential equation 

= s(l - u)N - fuN 0 < s < 1, 0 < / < 1 

Since N is constant then 

(*) u = — = s(l —u)—fu 



or 



du _ 

The Equilibrium requires that — = (J, or 

dt 



i.e., 



where 



il = S — (s + f)u 

s ~ (s + f)u* = 0 
* _ _£_ S J]_ 

U ~ s+f~ l+(s//) 

— = ^^>0 ^1 = ^^<0 



9s (s+/y ds (s+fy 

In other words, the equilibrium unemployment rate- the natural rate in this model- rises as more 
individuals enter the unemployment pool to actively search for a job and falls when the job-finding rate becomes 
greater. But this simple model says more than this about the equilibrium (natural) level of unemployment. It says 
that the level of u* occurs because individuals need to seek alternative unemployment and that the search for a 
new jobs takes time. 

The time path is readily found by solving the differential equation (*). If u(0) = u 0 , then 

u(t) = u* + (ii 0 - u*)e< s+ ^, u* = — 

s+f 

Since both s and / are positive, then this solution implies that unemployment tends to its equilibrium value over 
time. 

In this model, the focus is on the level of unemployment. Of course, if N is fixed, then the 
employment, E, is simply 

E = (1 - u)N 

or 

e = E/N = (1 - u), 

where e is the employment rate. In order to lay the foundation for other dynamic theories, it is worth noting that at 
any moment of time there will be an unemployment rate of u — U/N, and a vacancy rate of v — V/N. Since N 
is constant throughout, we can concentrate on the rates u, v and e. 

At any moment of time there will be an unemployment rate u and a vacancy rate v, where those 
individuals who are unemployed are attempting to match themselves with the available vacancies. Since we have 
assumed that the number of jobs is matched by the number of those seeking employment, then u — v., the 
problem is one of matching the unemployment to the vacancies. Accordingly, the literature refers to the matching 
rate or the exchange technology (Mortensen, 1990) [23]. In other words, the unemployment and the jobs that 
employers are seeking to fill are inputs into the meeting process. Let this be denoted m(u, v). 

Given m(u,v), then for such a meeting to take place must either be some unemployment or some 
vacancies. More formally m(0, v) — m(u,0) — 0. Furthermore, the marginal contribution of each input is 
positive,i.e., dm/ du > 0 and dm/ dv > 0. Following Diamond 1982 [14], it is further assumed that the 
average return to each input is diminishing, i.e., m/u and m/v diminishes with u and v, respectively. 
Finally, and purely for mathematical convenience, we assume that m(u, v) is homogeneous of degree k, so that 

m(u, v) = u k m{l,v/u) 

Using this analysis we can write the change in employment as the total match Nm(u, v) minus those losing a job 
s(l - u)N, i.e. 
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Jt = ^fP = Nm ( u > v ^ ~ s 0- ~~ U ^ N 

or 

— — e — m(u, v) — se 

dt V > J 

Although the time path of employment, e(t), must mirror the time path of the unemployment rate , u(t), since 
e — 1 — u, the present formulation directs attention to the matching rate m(u, v). 

In general (Mortensen 1990, [23]), the equilibrium hiring frequency, m(u,v)/u, is a function of the 
present value of employment per worker to the firm ,q, and the employment rate , e. This can be established by 
noting that 

, N m(u,v) u k m(l,v/u) -„ , N nJ/-1 / % i ^ % 

(**) — ^ = = it m(l, v/u) = (1 - e) fc 1 m(l, v/u) = e). 

The hiring function h(<7, e) is a function of q since the value of v/u in (**) is determined in equilibrium. In 
equilibrium, the return on filling a vacancy (mq/v) is equal to the cost of filling vacancy, c, i.e., 

771 (U, l/) 

[-^-Ar = c 
v 

that is 

(l-e) k - 1 q= cv/u 

^ CJ q m(l,v/u) 

so, the hiring frequency is related to both q and e. Furthermore, we can establish from this last result that h q > 0 
and ft e < 0 if k > 1 and ft e > 0 if k < 1. Hence 

m(u,v) , , . , _ 

u = h(q, e) h q > 0, 

System 1 

(h e < 0 if k> 1 
U e > 0 i/ < 1. 



m(it, = u/i(q, e) = (1 — e)/i(q, e) 

which in turn leads to the following equilibrium adjustment equation 

e = (1 — e)/i(q, e) — se 

The profit to the firm of hiring an additional worker is related to q and the employment rate e, e.i, n(q,e), and 
it will be different for different models of the labour market. The profit arises from the difference in the marginal 
revenue product per worker, MRP L , less the paid wage, w. If we denote the MRP L by x(e), then n(q,e) — 
x(e) — w. However, the future profit stream per worker to the firm is 

rq = x(e) — w — s(q — k v ) + q 
where rq represents the opportunity interest in having a filled vacancy and k v is the capital value of a vacant job 
,i.e., the present value of employment to the firm is the profit from hiring the worker less the loss from someone 
becoming unemployed plus any capital gain. 

Since in equilibrium no vacancies exist, then k v = 0 and so 

rq = n(q,e) - sq + q 
q = (r + s)q -n(q,e) 

To summarise , we have two differential equations in e and q i.e., 
System 2 

e = (1 — e)/i(q, e) — se 
q = (r + s)q - n(q, e) 

Whether there is a unique equilibrium, it rests very much on the degree of homogeneity of the match 
function i.e., the value of k in m(u,v) — u fc m(l,-) and the productivity per worker x(e). 
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III. Bifurcation Phenomena in Economics Literature 

Over the last three decades, the theoretical research in macroeconomics has moved from comparative 
statics to dynamics, with many such dynamical models exhibiting nonlinear dynamics. Bifurcation analysis is a 
key tool for the analysis of dynamic systems in general and nonlinear systems in particular. It has been widely 
used in mathematics and engineering and although it is a relatively new in economics, the interest for this type of 
analysis has been increasing because it provides information about the occurrence and changes in stability of fixed 
points, limit cycles and other solution paths, it helps model these changes and transitions from stable to unstable 
case or vice versa as some parameters change. Most common types of bifurcations encountered in economic 
analysis include Saddle Node, Transcritical, Pitchfork, Flip and Hopf bifurcations (William A. Barnett, 2008). 
Each of these bifurcations produces a different type of qualitative change in dynamics. 

The first paper on Hopf bifurcation belongs to Poincare' (1892) and the first formulation of the 
theorem to Andronov (1929). Another important theorem on the existence of Hopf bifurcation appeared in Hopf 
(1942). Its existence has been proved in many economic models. For example, Torre (1977), observed the 
existence of a limit cycle associated with a Hopf bifurcation boundary in a study on Keynesian systems. Benhabib 
and Nishimura (1979) explored a multi-sector neoclassical optimal growth model and proved that a closed 
invariant curve might be a result of optimization. Hopf bifurcations were also found in other studies such as 
Aiyagari (1989), Benhabib and Day (1982), Benhabib and Rustichini (1991), Gale (1973). Recent findings belong 
to Barnett and He (1999, 2001, 2002, 2004, 2006, 2008), who found bifurcation boundaries in a Bergstrom 
continuous -time model of the UK economy and to Leeper and Sims with their Euler-equations model of the 
United States economy. 

Barnett and Duzhak (2008, 2010) analyzed the bifurcation phenomenon using a closed economy New 
Keynesian model and they found both Hopf and Period Doubling bifurcations. Despite the growing interest in 
bifurcation analysis of economic systems, the literature on this subject is still immature and needs an extensive 
study for a comprehensive understanding. The present paper studies the existence of Hopf bifurcation for the 
model described by Ronald Shone (2002), by considering different settings for the matching function denoted 
m(u, v). 

A Preliminary Review of Hopf Bifurcation Theory 

According to Yuri Kuznetsov ( 2004), if we consider the following system of two differential equations 
depending on one parameters. 

Theorem 1 [20] Suppose a two dimensional system 

^ = f(x, a), x £ R 2 , a £ R 1 , (1.1) 
with smooth / , has for all sufficiently small \a\ the equilibrium x — 0 with eigenvalues 

^uO) = MO) + ico(oc), 

where ju(0) = 0, o)(0) = a> 0 > 0. 

Let the following conditions be satisfied: 
(B.l) Zi(0) & 0, where ^ is the first Lyapunov coefficient; 
(B.2) /i'(0) ±0. 

Then, there are invertible coordinate and parameter changes and a time reparameterization transforming 
(1.1) into 

^(pj = (l -p\ (pj ± ( yi 2 + y 2 2 ) (p + 0(||y|| 4 ). 

We can drop the 0(| |y| | 4 ) terms and finally arrive at the following general result. 
Theorem 2 ([20] Topological normal form for the Hopf bifurcation) 
Any generic two-dimensional, one-parameter system 

x=f{x, a), 

having at a — 0 the equilibrium x — 0 with eigenvalues 

h,2 («) = ± ia) o> w o > 0, 
is locally topologically equivalent near the origin to one of the following normal forms: 

£(£) = (i -J)(^)±(yi 2 +y 2 2 )(S). 

Consider the following system of two differential equations depending on one parameter: 
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System 3 



X-^ — OCX-^ X 2 "t" S ' "t" ^2^) 

X2 "I - CK%2 ^ ' -^2(^1^ ^2^)' 

This system has three equilibriums x 1 — x 2 — 0 for all a with the Jacobian matrix 

(a — 1\ 



having eigenvalues A 1>2 = a + i. 



Figure 1 .Supercritical Hopf bifurcation 
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Supercritical Hof Bifurcation 



Source: authors' estimations 



For s = — 1, System 3 always has an equilibrium point at the origin. This equilibrium is a stable focus 
for a < 0 and an unstable focus for a > 0. At the critical parameter value a — 0 the equilibrium is nonlinearly 
stable and topologically equivalent to the focus. Some times it is called a weakly attracting focus. This equilibrium 
is surrounded for a > 0 by an isolated closed orbit (limit cycle) that is unique and stable. 

The cycle is a circle of radius p 0 (a) = yfa. All orbits starting outside or inside the cycle except for the 
origin tend to the cycle as t -» +00. This is an Supercritical Andronov-Hopf bifurcation, this being represented in 
Figure 1. 



Figure 2 . Subcritical Hopf bifurcation 





for a=0 

Subcritical Hopf Bifurcation 



Source: authors' estimations 



For s — 1 in system 3 we have a Subcritical Hopf Bifurcation ( presented in Figure 2). 

The system undergoes the Andronov-Hopf bifurcation at a = 0. Contrary to System 3 for s = —1, 
there is an unstable limit cycle in 3, which disappears when a crosses zero from negative to positive values. The 
equilibrium at the origin has the same stability for a ^ 0 as in System 3 for s — —1: It is stable for a < 0 and 
unstable for a > 0. Its stability at the critical parameter value is opposite to that in 3: It is (nonlinearly) unstable 
at a = 0. 



Remark 1 We have seen that there are two types of Andronov-Hopf bifurcation. The bifurcation in System 3, for 
s = — 1, is often called supercritical because the cycle exists for positive value of parameter a ("after" the 
bifurcation). 

The bifurcation in System 3, for s — 1 is called subcritical since the cycle is present "before" the 
bifurcation. It is clear that this terminology is somehow misleading since "after" and "before" depend on the 
chosen direction of parameter variation. 
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IV. Examples of Dynamical Systems and their Stability Analysis 

This section develops and analyses two dynamical systems resulting from using, in System 2, the 
following expressions for the matching rate. 

m 1 (u,v) = a • o m u a v a {u a + v a ) 
m 2 (u,v) = o m u a v x ~ a 



where a is the elasticity of the matches with respect to unemployment but also the elasticity of the vacancy filling 
rate with respect to the labor market tightness , a m represents the efficiency of the matching process and a is a 
scale parameter. The first setting for the matching rate when m 1 (u, v) = uv(u + v) was proposed and analysed 
by L. C. Holdon (2011). We extend that function m 1 (u, v) to a more general one, where a 6 (0,1] is a scale 
parameter which can be used to particularize every economy, and we keep the linear function n(q, e): = q ■ e. 

The second setting is widely used in economics to help model the inefficiencies that are seen in many 
markets where two agents seek out each other in order to come to an agreement (search frictions). Though, its 
main application is within the labour market. 

The matching function theory has helped to improve models related to business cycles and has become 
the most significant and important tool in analyzing the labour market, particulary with regards to macroeconomy. 
It is used by a variety of governments when trying to decide upon unemployment policy, and macroeconomic 
policy as a whole , its simplified implication being that the lower the level of search friction, the more efficient the 
relevant market is, and as such the lower the cost that is associated with the relevant markets pairings. 
The greatest succes and research in this sense belongs to Diamond, Mortensen and Pissarides, whose model is the 
most prominent macroeconomic tool which shows the relationship between the rate at which the unemployed are 
hired with the number of people looking for jobs, and the volume of vacancies available. 
Further on, the study will explore the ode systems from a mathematical view point. 



From introduction we know that U = 1 — e and denote q : = — . 

e: = xG(0,l) 
q: = y ER 



Namely, 
and by equation (**) we have 



miu.v) = a-o m {l -x) 3a y a (l + y a ) 
\m 2 (u,v) = <j m (l -x) -y 1-0 '. 



Following System 2, we obtain the o.d.e systems of: 
System 4 



fx = a ■ o- m (l - x) 3a y a (1 + y a )) - s ■ x -.= f x {x,y) 



fx = a ■ er m (l — x) 
\y = (r + s~) ■ y — x 



y :=/ 2 (x,y) 



System 5 

x = o- m (l -x) -y 1 -* 7 -s -x ■■= s-H^y) 

y = (r + S ).y- X .y :=g 2( xy > ) 

where (r, s) are parameters, with s £ (0,1) and rEl; 

x is the employment rate and y represents the labour market tightness. 

The Study of System 4 

Applying the continuous transformation (x, y 17 ) — > (*<y) we obtain for the first equation of System 4: 

x = a • o m (l - x) 3t7 y(l + y) - s • x, 

and for the second equation we obtain: 

y = a • y a ~ 1 • y = a • y a ~ 1 • [(r + s) • y - x • y] = a • [(r + s) • y a - x ■ y a ] = 
<r[(r + s) • y — x • y]. 

Therefore, the equivalent ode system with System 4 is: 
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System 6 

(x = a • <7 m (l - x) 3t7 y(l + y)-s -x := f x {x,y~) 
[j = a [(r + s) • y - x ■ y] ■= f 2 (x,y) 

For System 6, we have three equilibrium points: 



E 0 = (0,0), £i = (r + s,-|-|- ll+ 4 ;f +S \ 3 7) 

2 2 -\J a-tr m (l-r-s)- 50 ^ 



and 



1 1 4s(r + s) 
£ 2 = (r + s,- - + -• 1+ ^ ^j). 



2 2 J a-o m {l-r-sy 

The Jacobian matrix associated to system 4 is: 

f-3-a • o • o m • (1 x)^" 1 ) • Kl + - s a • a m • (1 - x) 3a (1 + 2y") N 
^ : = I — <ry cr[— x + (r + s)] 

where £" = (x, y) is an equilibrium point. 
Since S G (0,1), we have: 

Proposition 1 The equilibrium point E 0 is: 
(t) an unstable saddle point for r + 5 > 0. 
(it) an attractive stable node for r + s < 0. 

(tit) a nonlinearly attractive equilibrium by type saddle-node for r + 5 = 0. 
The Characteristic Equation associated to point E 0 is: 

A 2 - (-s + a ■ (r + s))X + (-s) • a • (r + s) = 0, 

with 

A x = (s + er(r + s)) 2 => /I4 = a • (r + s) and 1 2 = — s > 

because 0 < S < 1 => 1 2 < 0 for a11 s - 
For the equilibrium E 1 we obtain the following: 

rrt , o S(r+S) 

Tr/i = -3 • a • — — - - s 

l-r-s 

and 



with its characteristic equation: 

A 2_ ( _ 3 . f7 .^±£)_ s)A+( l. a . f7 . (1 _ r _ s) 3a /l+ 4 f +S \ 3ff (1 + 

Ii + 4s(r+s) , )) = 0. 

Denoting: 



p-.= 



1 V a!ll-r- s)* = ^ ■ a ■ a ■ a„ ■ (1 - r - . p. (1 + W . 
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Figure 3. Parameter portrait for System 6 for the equilibrium points 











in 






V 



TrJ_1=0 DetJ_1=0 

Source: authors' figure 

Proposition 2 The equilibrium point E { (see Figure3) is: 

(i) For zone I, it is a stable attractive node when \ 2 €E R , and stable attractive focus when \ 2 £ R, with 

Re\ 2 <Q. 

(ii) For zone II, it is a nonlinearly point of Hopf type. 

(iii) For zone III, it is an unstable repulsive node when /lj 2 £ R + , and an unstable repulsive focus when 

\ 2 £ R, with ReA { 2 > 0. 

(iv) For zone IV, it is a nonlinearly equilibrium by type saddle-node. 

(v) For zone V, it is an unstable saddle. 

From mathematical view point, the next result is the most important, it describes the existence of Hopf 
bifurcation in System 6, with serious repercussions for a country's economy. 

1 

Theorem 3 If V + S H = 0 and Detji > 0, corresponds to zone II. Then the equilibrium point 

3-c—l 

E l is a nonlinearly equilibrium point of Hopf type. 
Proof. We must verify the two conditions from Theorem 1 . 

First, we must calculate li (r + S + = 0^, which is the first Lyapunov coefficient, when V + S + 

1 1 
= 0. For obtaining U (r + S H = 0), we will use the formula from [181: 

h(X, y)| . '•= —UL + fly + fly + fly\ +— + f]y)- f%Ul + flfl + flyfyyl 

(r+i+_!_=0) lO 0) Q lO 

3-ct-I 

where 
System 7 

fx = a • er m (l — x) 3t7 y(l + y) — s • x := / 1 (x,y) 
|y = <7[( r + s) -y-x-y] :=/ 2 (x,y) 



and 



xxx : ~ ~~ aT 3 - 



Since O m , O, S > 0 and r + S + _ = 0=>r + S=£l=>/?>0 

co 0 = 2 • Ji • a • o- m • a • (1 - r - s) 3ff • /? • (1 + /?) > 0 



Since = /yy = 0 we have 



1 1 1 1 1 1 1 1 

^lC^»y)| i "T7 7 \fxxx ~l~./scyy]~l fxy (fxx fyy)> 

('•+s+3^3T=o) 16 (x)q 16 ' ^ 



hence 
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h(x,y)\ i 

(r+s-f 



= 0) 



(-1) • • a • [3 • a m ■ a • (3<r - 1) • (3a - 2)(1 - 

i i 



x) 3 ^" 1 )y(l + y) + 6-a-a m •(l-x)^ 3ff - 1 )]+ — — 3 • a 2 ■ o ■ o m ■ (1 - x)( 3ff-1 >(l + 

60q 16 

2y) • [3a • o m ■ (3a - 1)(1 - x)^" 2 )y(l + y) + 2 • a m • (1 - x) 3 *]}. 

l 



We recall that X = T + 5 = 



1-3-CT 

Zi(r + s + 



1 and <T m , (7 > 0, so, we deduce that 
1 

= 0) * 0 , /or 0< s < 1. 



3 • a- 1 



And the second condition from Theorem 1, since r + S H 

3-tr-l 

a(r+s+r-^—) d(r+s+— !—) 



0 and S £ (0,1):: 



5s 



9r 



Following Theorem 1 we deduce the existence of Hopf bifurcation. 

If inside of Theorem 3 we consider (X = o~ = o~ m = 1 then we obtain the next result, presented in 

1 

Holdon(201 1) (see [22]), with a clear picture which illustrates the Hopf bifurcation for V + 5 + - = 0, 
associated to the system: 

System 8 

fx = (1 — x) 3 • y • (1 + y) — s • x ■= / 1 (3f,y) 
iy = (r + s) • y — x • y := f 2 (x,y) 



Corollary 17fr + S + ~ = 0 awe/ DetJi > 0, fin System 8), corresponds to zone II. Then the equilibrium 
point E l is a nonlinearly equilibrium point of Hopf type, see Figure 4, and the Hopf bifurcation is subcritical. 

Figure 4: r + s + ~ = 0, see [22] 

Initial data : 
x=- 0.49099 
y=- 0.02536 

Browser: 
200 
200.01 
200.02 

For integration 
parameters 
we chose the 
Plick method : 
RK4. 




The value of the parameters: ( r=-1 ; s=0.5 ) 

Integragration parameters : Delta t=0.0005 ; T start=0 ; bounds=50 ; Max step=1 . 



Source: authors' figure 

For the equilibriums E 2 we obtain the following: 



s(r + s) 
TrJ 2 = — 3 • a • s 



1 — r — s 



and 



DetJ 2 = - - • a • a • a m • (1 - r - s) 



3a 



L 4s(r+s) _ _ L 4s(r+s) . 

a] a-a m (l-r-sf a ^ a] a-a m (l-r-sf a) ' 
s(r+s) 



with its characteristic equation: A 2 — (— 3 • O • — s)X + (— — • Ct • a ■ O m ■ (1 
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s) 3a L 4s(r+s)~ I 4s(r+s) ~ = 

Aj a-tr m (l-r-s) 3ff v aJ a-ff m (l-r-s) 3ffyy 

Figure 5. Parameter portrait of system 6 for the equilibrium point E 2 
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Source: authors' figure 

Proposition 3 The equilibrium point E 2 (see Figure 5) is: 

(i) For zone I, it is an unstable saddle. 

(ii) For zone II, it has the eigenvalues 



2 \ a-cx m (\-r-s) \ a-a m (l-r-s) 

because T + S H = 0 and it is an unstable saddle. 

3-cr-l 

(iii) For zone III, it is an unstable saddle. 

(iv) For zone IV, it is a nonlinearly equilibrium by type saddle-node. 

(v) For zone V, it is is an unstable repulsive node when A { 2 £ R + , and an unstable repulsive focus when 

\ 2 £ R, with Re\ 2 > 0. 
The Study of System 5 

For System 5, we have two equilibrium points: 

l 

E 2i0 = (0,0), and E 2fl = (r + s, 5(r ^ 

The Jacobian matrix associated to System 5 is: 

f-W^-s a m (l-x)(l-a)y"- 
/: = I — y r + s — x 

where E = (x, y) is an equilibrium point. 
The study of the equilibrium point E 2 Q : 

Since S G (0,1) we have: 
Proposition 4 The equilibrium point E 2 0 is: 

(i) an unstable saddle point for r + s > 0. 
(z7) an attractive stable node for r + s < 0. 

(iii) a nonlinearly attractive equilibrium by type saddle-node for r + s = 0 
The Characteristic Equation associated to point E 2 0 is: 

A 2 -rX + (-s)(r + s) = 0, 
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with 

Aa = (r + 2s) 2 => X x = r + s and A 2 = —s, 

because 0 < S < 1 => A 2 < 0 for all S. 
For the equilibrium E 2 l we obtain the following: 

TrJ 2il = 



and 

with its characteristic equation: 

A 2 + 



(1 — r — s) 
DetJ 2:1 = (1 - a)s(r + s). 

-A + (1 - a)s(r + s) = 0. 



(1-r-s) 

Corollary 2 There isn't Hopf bifurcation for System 5, since TvJ 2 i ^ 0 for any S G (0,1). 

Figure 6. Parameter portrait for the system 5 for the equilibrium point £ 2 ,i 




0; 

DetJ_2,1=0 

Source: authors' figure 

Propositions The equilibrium point E 2 l (see Figure 6) is: 

(i) For zone I, it is an unstable repulsive node when A\ 2 £ R + , and an unstable repulsive focus when A\ 2 £ R, 
with /^g/lj 2 > 0. 

(ii) The second zone doesn't exists, since TrJ 2l ^Q for any 5e(0,l). 

(iii) For zone III, it is a stable attractive node when A\ 2 e R , and a stable attractive focus when ijj^R, 
with /te/lj 2 < 0. 

(iv) For zone IV, it has the eigenvalues A\ = 0 and A\ < 0 and it is a nonlinearly attractive equilibrium by 
type saddle-node. 

(v) For zone V, it is an unstable saddle. 

V. Conclusions 

The purpose of this study is to develop dynamic models in order to point out a set of facts regarding job 
flows and unemployment. The main results of this research are the theorems on the existence of Hopf bifurcation. 
The proposed models extend Ronald Shone's approach to the dynamics of unemployment by providing two 
different settings for the matching rate m(u,v). One of the settings is the matching function proposed by the 
distinguished Nobel Laureates Diamond, Mortensen and Pissarides, the other one being a personal attempt to find 
another function that could provide stability to a dynamic system regarding unemployment. The stability of the 
resulting systems is analyzed by using the software package of xpp.exe. 

The most important findings are: the ode system resulting from the function that we have proposed has 
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three equilibrium points among which the second zone of the equilibrium point E [ is a nonlinearly point of Hopf 
type describing the existence of a Hopf bifurcation; the ode system resulting from using the matching function 
proposed by Diamond, Mortensen and Pissarides has two equilibrium points among which a nonlinearly attractive 
equilibrium of saddle-node type. The system resulting from using the setting proposed by Diamond, Mortensen 
and Pissarides is more appropriate to the economic reality, this study being another proof that their function is the 
most significant and important tool in analyzing the labour market dynamics. But from mathematical point of 
view, the system that we proposed is much more interesting because it allows the existence of Hopf bifurcation. 
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